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We study the inverse Compton scattering of the CMB photons off nonthermal high-energy elec- 
trons. In the previous study, assuming the power-law distribution for electrons, we derived the 
analytic expression for the spectral intensity function in the Thomson approximation, which 
was applicable up to the photon energies of a; < O(GeV). In the present paper, we extend the 
previous work to higher photon energies of a; > O(GeV) by taking into account the terms dropped 
in the Thomson approximation, i.e., the Klein-Nishina formula. The analytic expression for /(oi) 
is derived with the Klein-Nishina formula. It is shown that has a "knee" structure at uu = 

O(PeV). The knee, if exists, should be accessible with gamma-ray observatories such as Fermi-LAT. 
We propose simple analytical formulae for which are applicable to wide photon energies from 
Thomson region to extreme Klein-Nishina region. 
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Keywords: cosmology: cosmic microwave background — cosmology: theory — galaxies: clusters: general — 
radiation mechanisms: thermal — relativity 



I. INTRODUCTION 

The inverse Compton scattering is one of the most fun- 
damental reactions which have variety of apphcations to 
astrophysics and cosmology. They are, for example, the 
Sunyaev-Zeldovich (SZ) effects T] for clusters of galax- 
ies (CG), cosmic-ray emission from radio galaxies and 
clusters of galaxies Q, and radio to gamma-ray emission 
from supernova remnants [1, Therefore, theoretical 
studies on the inverse Compton scattering have been 
done quite extensively for the last forty years, start- 
ing from the works by Jones [6l|, and Blumenthal and 
Gould to the recent works, for example, by Fargion[8i], 
Colafrancescofgj [l0| . and Pctruk^ll]. 

In particular, remarkable progress has been made in 
theoretical studies for the SZ effects for CG. Wright 
and Rephaeli[l^ calculated the photon frequency redis- 
tribution function in the electron rest frame, which is 
called as the radiative transfer method. On the other 



hand, Challinor and Lasenby[l4| and Itoh, Kohyama, and 
Nozawa[15j solved the relativistically covariant Boltz- 
mann coUisional equation for the photon distribution 
function, which is called the covariant formalism. Very 
recently, however, Nozawa and Kohvama|,16l] showed that 
the two formalisms were mathematically equivalent in the 
approximation of the Thomson limit. 

On the other hand, Nozawa, Kohyama and Itoh(l7l.[l8j 
(denoted NKI hereafter) extended the formalism ob- 
tained by T&l to the case of high-energy electrons in the 
Thomson approximation. Assuming the power-law dis- 
tribution for electrons, the analytic expression was de- 
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rived for the spectral intensity function I{uj), where w is 
the photon energy. This extension was particularly in- 
teresting for the analysis of X-ray and gamma-ray emis- 
sions, for example, from radio galaxies [2| and supernova 
remnants 0, where the inverse Compton scattering of 
the CMB photons off nonthermal high-energy electrons 
plays an essential role. 

In the present paper, we extend the formalism obtained 
by NKI to higher photon energies beyond the Thom- 
son approximation by calculating the Klein-Nishina for- 
mula. The present extension will be particularly im- 
portant to study the inverse Compton scattering of the 
CMB photons for the photon energies ui > O(GeV). We 
will derive the analytical expression of /(w) for the case 
of the power-law electron distribution with the Klein- 
Nishina formula. It will be demonstrated that /(w) has 
a "knee" structure a.t w — O(PeV). The knee, if exists, 
should be accessible with gamma-ray observatories such 
as Fermi Gamma-raySpace Telescope's Large Area Tele- 
scope (Fermi-LAT) jl9j . We propose simple analytical 
formulae for /(w) which are applicable to the entire pho- 
ton energies. 

The present paper is organized as follows: In Sec. II, 
we summarize the results of the inverse Compton scat- 
tering in the Thomson approximation, which was ob- 
tained by NKI. We calculate the rate equation and show 
the analytic expression for the spectral intensity function 
dI{uj)/dT in the Thomson approximation. In Sec. Ill, we 
extend the formalism of Sec. II with the Klein-Nishina 
formula, and derive the analytical formula which is ap- 
plicable to higher photon energies of > O(GeV). In 
Sec. IV, we discuss various other cases for the boundary 
parameters of the electron distribution function which 
are not studied in Sec. III. Finally, concluding remarks 
are given in Sec. V. 
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II. HIGH-ENERGY INVERSE COMPTON 
SCATTERING IN THOMSON APPROXIMATION 



terms of the power-law distribution function of three pa- 
rameters: 



A. Rate equations 

In the present section, we summarize the results ob- 
tained by NKI 17] in order to make the present paper 
more self-contained. The rate equation for the photon 
distribution n{ui) and the spectral intensity function 
were derived under the following Thomson approxima- 
tion: 



7— < 1, 



(1) 



where uj is the photon energy, 7 is the Lorentz factor, and 
m is the rest mass of the electron. Throughout this pa- 
per, we use the natural unit h = c = 1, unless otherwise 
stated explicitly. They are given as follows: 



dn{(jj) 
dr 



dsPi^TA{s) [n{e''xT) - n{xT)] , 



J —00 

POO 

Pi,ta{s) = / d-fpe{j) Pta{s,j) , 



'7* 

r = nea-rt , 

UJ 

XT 



MB 



(2) 
(3) 

(4) 

(5) 
(6) 



The subscript TA denotes that the Thomson approxi- 
mation is used in the calculation. In above equations, 
I{uj) — IoXj^n{uj), Iq — (kBTcMB)^ /2tt^ , is the elec- 
tron number density, ctt is the Thomson scattering cross 
section, and Pi,ta(s) is the probability distribution func- 
tion for the photon of a frequency shift s defined by 
e'* = x'^/xT. 

In the present paper, we are interested in the inverse 
Compton scattering of the CMB photons by high-energy 
electrons, where 7 ^ 1 is assumed. In this approxima- 
tion, the minimum value 7* of the 7-integration is ap- 
proximated by 



7^ 



(7) 



Then, we havefTTj the analytic expression for Pta(si7) 
as follows: 



Pta(s,7) 



/2 



^2 



(8) 



In Eq. Pe(7) is the electron distribution function, 
which is normalized by J d'ype{'^) — 1. High-energy elec- 
trons in the supernova remnants and active galactic nu- 
clei, for example, are most likely nonthermal. Therefore, 
it is standard to describe the nonthermal distribution in 



Veil) 



N^l " , 7min < 7 < 7max 

, elsewhere 



(9) 



where Nj is the normalization constant, a is the power 
index parameter, 7niin and 7max are parameters of mini- 
mum and maximum values for 7, respectively. 



B. Spectral intensity function 

In NKI, the rate equations of Eqs. ([2]) and ([3]) were 
solved for the Planckian initial CMB photons with the 
power-law electron distribution of Eq. For the inverse 
Compton scattering of the CMB photons by high-energy 
electrons, we are interested in high-energy spectrum such 
as X-rays (~ keV) and gamma-rays MeV). Therefore, 
one can safely assume a;T ^ 1 for scattered photons, 
where the second terms in Eqs. ([2]) and ([3]) can be ne- 
glected. For the 7-parameters, we assume the following 
condition: 



1 < 7min < 7 



(10) 



Under these assumptions, Eqs. ^ and ^ were solved 
analytically. One obtained as follows: 



d/TA(t^) 



dr 



3Xt 



dt- 



e* - 1 (T + 1 



cr + 3 Vcr + 3 t / t 



+ 



6(^2 



4(7+11) 



(CT-Hl)(a + 3)2(a + 5) 



-(<T-l)/2 



Xt ^((J+3)/2 

dt— 



XT 
2 

min 



47 



(11) 

(12) 



Note that — (c— l)7^ijj^ was used under the condition 



of Eq. (fTO|) . 

For Xt ^ 1, Eq. (fTTT) was further simplified as follows: 



dlTAji^) _ 6((7- l)(g2 +4(7 + 11) 
dT " ° ((7 + l)(a + 3)2(a + 5) 



C(^) 



r(2) Jo 



dt- 



1 ' 



(13) 
(14) 



where C,{z) is the Riemann's zeta function. In Fig. 3 of 
NKI, dlTp,{oj) / dr of Eq. (ITT]) were plotted as a function 
of Xt for typical cr-values. It was shown that Eq. 
was a good approximation of Eq. pT|) for > 10. 
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III. HIGH-ENERGY INVERSE COMPTON 
SCATTERING WITH KLEIN-NISHINA 
FORMULA 

A. Rate equations 

In the present section, we extend the formaUsm de- 
veloped in NKl[l3] to higher photon energies where the 
Thomson approximation breaks down. In order to ac- 
comphsh the task, one needs to relax the condition of 
Eq. ([1]). Instead, one needs to keep the full (so called 
Klein-Nishina) amplitude for the Compton scattering 
process. 

Very recently, Poutanen and Vurm 20] made an exten- 
sive study on the Compton scattering process with the 
general Klein-Nishina cross section. The formalism can 
be applicable to higher energies. On the other hand, it 
was shown by Nozawa, Kohyama and Itoh^lj that their 
formalism was equivalent to the NKI formalism in the 
Thomson approximation. In the present paper, there- 
fore, we start with the formalism of [l^l and derive the 
analytic expression for the spectral intensity function. 

The source function is defined by Eq. (174) of [l^ and 
the rate equations for the photon number density and the 
spectral intensity function are given as follows: 



dn{ijj) 

dr 
dljuj) 

dr 



rfsPi^KN(s) ["(g^xt) - n(a;T)] , (15) 
dsPi,KN(s) [I{e-'XT) - /(xt)] , (16) 

(17) 

(18) 
(19) 



^'i,kn(s) = / d7Pe(7)PKN(s,7) , 

^'kn(s, 7) = 16^^^^^ y ' 



X 



The subscript KN denotes that the Klein-Nishina cross 
section is used in the calculation. In Eq. pT)) . 7* is the 
minimum value of 7. In Eq. (|18l) . /imin and /J.max are 
minimum and maximum values of the cosine of the scat- 
tering angle, and the function i?o is given by Eq. (138) 
of [20]. 

Let us consider the high-energy (x » 1) scattering of 
the CMB photons {x' ^ 1) by extremely high-energy 
electrons (7^1, therefore /3 « 1). It should be empha- 
sized here that all conditions made in Sec. II are retained 
except for Eq. ([Ij. In this approximation, one finds 



7* 



Mmin — 1 5 

M 



a 



27(7 - x) 



(20) 

(21) 
(22) 

(23) 



Then, Eq. p8)) can be approximated as follows: 
3 



fkN(s,7) = 



I672 



2e''h+xe'' {hil - x) - hij)} 



-— {/3(7-x)-/3(7)} 

X 

-l{I,{j-x)-h{j)} 

+ \{il-x)l5{-f-x) -7/5(7)} 



x^ 

X 



{I,{^-x)-Ieh)} 



(24) 



where Ii, l2{(,), ■ ■ ■ , Ie{^) are the integrals defined by 



h = 



dji- 



_i '(l + e2''-2e''^)i/2 ' 
1 

1 
1 

_1 '"(1-^)2(^2+^)3/2' 

rM„.a. 1 
> 1 ''(l-M)(e2+r)3/2' 



(25) 
(26) 
(27) 
(28) 
(29) 
(30) 
(31) 



Their elementary integrals can be performed. One has as 
follows: 



2s -1 1/2 



1 + e"- {1-2(1 -a)e" + e2^} 



(32) 



/2(e) 



^2-1 V2 



(^2 _ 1)3/2 



'1 + 1/^ 



(33) 



htO = - 



. (34) 



(35) 



4 



ve 



where 



(36) 



(37) 



(38) 



Under assumptions of high-energy scattering, namely, 
s < -1, > 1, and a = 0(1), Eqs. (I221)-(I371) are further 
simphfied by keeping only leading-order terms. One has 



/i 2 - a + 0{e'') , 

^(i-f).o(r3), 



/2(e) 



i In I + o(r^). 



o(r 



(39) 
(40) 

(41) 

(42) 

(43) 
(44) 



Inserting Eqs. (139])-(|44]) into Eq. ([24]), and perform- 
ing a lengthy but straightforward calculation, one finally 
obtains the following analytic expression for Pkn(s,7): 



1/1 1 
/5(0 = ^(l-f)+0(r^), 



Pkn(s,7) = 



3e3^ 
3272 

2 



1 



7^(7 — xy 



o3|s|/2 



J«l/2 + 8eH^I/2-|-^^ln^ 
7(7 — x) 7(7 — x) 2 



4x2 



7(7 - 2^) 



-|s|/2 



(45) 



where a is given by Eq. (P^. It should be remarked 
that Eq. P5|) returns to Eq. ([5]) by taking the Thomson 
limit J ^ X and a — > e' "1/272. In deriving Eq. (|45p . the 
formula 

Pkn(s,7) = e''^KN(-s,7) (46) 
was also used in order to derive the expression for s > 0. 

B. Spectral intensity function for x > 7min region 

In the present subsection, we now calculate the spectral 
intensity function with Eq. (HSl) . In order to proceed the 
calculation, we use Pe{l) of Eqs. ^ and (fTO)) with an 
additional condition x > 7min- For simplicity, we also 
put 7max -> 00. (We study x < •j.nin case and finite 
case in Sees. IV-A and IV-B, respectively.) 



Let us first introduce new variables u, t and X as fol- 
lows: 



7 = x{l + u) . 
t — e^XT , 
X = Ocmbx , 



(47) 
(48) 
(49) 



where ^cmb = k^TcMB / itt-c^ ■ Then, we rewrite Eq. (^5 
with these variables. One obtains as follows: 



3^CMB 



1 



1 



32x4(1 -I- u)2 

2 t 
^Xu{l + u) 

4 t 

4^2 1 



X2u2(l+u)2 

{\n4:Xt + \nu{l + u)} 
t 



u{l + u) Xm2(1 + u)2 



(50) 



Inserting Eqs. ([9]) and (f50| into Ec^. ([TTl) . and perform- 
ing a straightforward calculation, one finally obtains as 
follows: 

^■^kn(^) 3/0, 

J = -;3r(CT - Ij (7minC'CMBj X QKN{X,a), 

cLt 32 

(51) 



)kn(^,<t)= / dt- 



Q{X,t,u) = 



1 



duQ{X,t,u), (52) 



1 ln4Xi 



(l + u)'^+2 ' ' 2Xt Xt 
1 4 ln[u(l + u)] 

1 A 1 



Xt\ xt J u2(l + u)'^+4 




(53) 



(54) 



where is the minimum value of the u-integration. In 
deriving Eq. (1521) . the Planckian distribution function 
was used for the initial CMB photons, and the second 
term in Eq. ([T6| were neglected, because xt ^ 1 is valid. 

Equation (|52p is one of the main results of the present 
paper. The spectral intensity function /kn(w) is ex- 
pressed by the double integral form. Therefore, the 
computation can be done extremely fast compared with 
the full numerical calculation of multi-dimensional inte- 
grals which involves numerical cancellations among vari- 
ous terms. 
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C. Approximate forms for spectral intensity 
function 

The expression for dlKyi{^)/dT can be further simph- 
fied for X <C 1 and X ^ 1 regions. For example, can 
be approximated by 



1 



2VXt 
1 



for X < 1 



for X > 1 



(55) 



4Xt 

where t = 0(1) is assumed. 

(i) X <C 1 region 

First, for X 1 region, the u-integration can be done 
by expanding Eq. (j53p in powers of l/u. Then, the t- 
integration can be also done, and one obtains the follow- 
ing analytic expression: 



dT 



Iq Xrj, 



-(<T-l)/2 



6{a -l){a^ +Aa + 11)^ fa + 5 

r 



(a + l)((7 + 3)2(a + 5)' V 2 
6(cr - l)(cr2 + 6cr + 16)^ i^cr + 6 
(a + 4)2(a + 6) 



2 

(7 + 6 



-o{x) 



(56) 



It is needless to mention that the condition ^ 1 is 
satisfied in this region. The first term in Eq. (1561) coin- 
cides with Eq. (jl3p , and the second term is the first-order 
correction to the Thomson approximation. It should be 
remarked that Eq. ([56)) was originally derived by Blu- 
menthal and Gould fTil in their Eqs. (2.76) and (2.77). 

(ii) X ^ 1 region 

On the other hand, for X 1 region (so called extreme 
Klein- Nishina limit), the expression for diKN {ijj) / dr can 
be also simplified. One obtains 

= -^^(cr - Ij (7min6'cMBj ^ Q^j^{X,(T), 

(57) 

Q'i^{X,a) = 4C(2) [\nX + C(a) + 1 - - C,] , (58) 

where C£;=0.5772 is the Euler's constant, and Ce ~ 
0.5770. We show the explicit derivation of Eq. (1551) in 
Appendix A. Note that the original expression was de- 
rived by Blumenthal and Gould?| in their Eq. (2.88). 
One can rewrite Eq. (2.85) in [7] as 



Cbg(^t) 



du 



V+3 



1-K21n2. 



(59) 



On the other hand, according to Eq. (jAlOp we have 

1 1 



C{a) 



du 
u 



{l + u) 



- 1 



(T + 3 



(1 + It)" 

71—1 ^ 

21n2-^- 



(60) 



It is clear that Cbg(o') corresponds to n = 1 of Eq. ([^ . 
In the present calculation, we choose the integer n w cr+S 
in order to minimize the integration in Eq. (1601) . 

It is important to note that Eq. (|60p is further simpli- 
fied for a = integer cases. One has 



2 "^'^ 1 

C(a) = — --l + 21n2-5]-. 



(61) 



because the integrand in Eq. (ISO)) vanishes by choosing 
n = cr -|- 3. The analytic expression is also obtained for a 
= half-integer cases as follows: 



C{a) 



1 



1 + 4 In 2 



^ 2i 

i=l 



1 



1 ' 



(62) 



where [a] denotes the integer part of a. The explicit 
derivation is shown in Appendix B. Similarly, it is shown 
that C((t) is written in analytic forms for a = arbitrary 
fractional numbers. In Appendix B, we show their ex- 
plicit forms. These analytic forms will be extremely con- 
venient to use compared with the original integral form. 



D. Knee structure of spectral intensity function 

In the present paper, we have derived the expres- 
sion for the spectral intensity function dI{Lo)/dT in the 
Klein- Nishina region, where the Thomson approximation 
breaks down. Equations pip and ([5T|) will provide reli- 
able theoretical data of the high-energy inverse Comp- 
ton scattering from Thomson region to extreme Klein- 
Nishina region. It is known from Eqs. (|56|) and (|57|) that 
the slope of the spectrum changes as 



dl[uj) ' 



^"-^y^ for X < 1 
for X > 1 



(63) 



The transition of the slope is expected to occur dX X 1, 
which corresponds to 



10^5 eV . 



(64) 



However, the transition can be checked only through per- 
forming the full calculation with Eq. ([55)) . 

In Fig. 1, we have plotted dI{u!)/dT as a function of 
X for a typical value a = 2.5. The solid curve is the full 
calculation of Eq. (I5T)) with Eq. ()52l) . The dashed curve 




FIG. 1: Plotting of dI{u))/dT as a function of X for o = 
2.5. The solid curve is the full calculation of Eq. (jSip with 
Eq. ((52} . The dashed curve and dash-dotted curve correspond 
to the approximations of Eqs. H56p and (I57p . respectively. The 
dotted curve is the Thomson approximation of Eq. (|13|l . 



and dash-dotted curve correspond to the approximations 
of Eqs. and ([57)) . respectively. The dotted curve is 
the Thomson approximation of Eq. p^ . The full calcu- 
lation clearly shows the knee structure at X « 1, which 
corresponds to a; w 10^^ eV. The knee, if exists, should be 
accessible with gamma-ray observatories such as Fermi- 
LATfli. 

On the other hand, the knee structure in the cosmic ray 
spectrum has been known since 1958 [l^. For the recent 
review, for example, see Bliimer, Engel and H6randel[23']. 
The main sources of the cosmic ray are charged particles 
such as protons. He and heavier nuclei. The observed 
energy spectrum shows a clear bending at around 10^^ 
eV. The energy spectrum follows a power law dN/dE cx 
E^°-. The power index increases from a=2.7 to a=3.1 at 
E « 10^^ eV, which gives the increment of Aa=0.4. The 
origin of the knee structure, however, is yet unsettled. 

In the case of the inverse Compton scattering of the 
CMB photons, the increment of the power index is given 
by Aa ^ {a + l)/2 1.75 for a = 2.5. Although the 
position of the knee is same, the change of the slope is 
quite large compared with the case of the cosmic ray 
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I J 

10^ 10*^ 10^ lO'^ lO'^ 10^'' 10^^ 

cj (eV) 



FIG. 2: Plotting of dl(uj)/dr as a function of uj for 7min = 10'^. 
The solid curve, dash-dotted curve and dashed curve corre- 
spond to a = 2.5, 3.5 and 4.5, respectively. In the calcula- 
tion, Eq. pip is used in the Thomson region {uj < 10^ eV), 
and Eq. H5ip with Eq. (|52p is used in the Klein-Nishina re- 
gion (a; > 10^ eV). The two curves coincide at the boundary 
region. 



spectrum. 

In Fig. 2, we have shown a spectral intensity function 
dI{uj)/dT for a wide energy range from uj = 10 eV to 10^^ 
eV. A typical value 7min = 10'^ is used for an illustrative 
purpose. The solid curve, dash-dotted curve and dashed 
curve correspond to cr = 2.5, 3.5 and 4.5, respectively. In 
the calculation, Eq. (jlip is used in the Thomson region 
(w < 10^ eV), and Eq. 1^ with Eq. ^ is used in 
the Klein-Nishina region (w > 10^ eV). The two curves 
coincide at the boundary region. It should be remarkable 
that two theoretical formulae have the prediction power 
in the energy range of the twenty orders of magnitude. 



IV. MISCELLANEOUS CASES WITH 
KLEIN-NISHINA FORMULA 

A. Spectral intensity function for x < 7min region 

In deriving the spectral intensity function with the 
Klein-Nishina formula in Sec. Ill, we have restricted our- 



7 



selves to the case where the photon energy x satisfies the 
following condition (denoted Case I hereafter): 

7inax oo • • • Case I . (65) 



In Case I, it is needless to mention that Eq. (|5ip was not 
applicable to the photon energy region x < 7min- For 
the X < 7niin region, however, Eq. ([FT]) was smoothly 
connected to Eq. (|13p of the Thomson approximation. 

In the present subsection, on the other hand, we derive 
the spectral intensity function with the Klein- Nishina for- 
mula for the X < 7inin region under the following condi- 
tion (denoted Case II hereafter): 



^ 7min < 7 max ; 7max — > oo • • • Case II . 



(66) 



In order to accomplish the task, we first rewrite the 
rate equation of Eq. (jl6p with the explicit lower (7rnin) 
and upper (7max) boundary parameters as follows: 



dljuj) 



/OO 

d7Pe(7)PKN(s,7) , (67) 



lnax(7^^. ,7min) 



-^inax — lll[4'7max('yniax ■2')] ■ 



(68) 



where 7* is given by Eq. (|20p , and we dropped the second 
term in Eq. It is clear that Eq. (|67l) returns to 

Eq. under the condition of Case I. 

Under the condition of Case II, Eq. (|67l) becomes 



dljcu) 
dr 



dsl{e "xt) / d7Pe(7) Pkn(s,7) 



+ / dsI{e~'xT) d7Pe(7)^KN(s,7), (69) 

Qin " Tmin 



Smin = ln[47min(7min " x)] 



(70) 



Repeating the same procedure done in Sec. Ill, one finally 
obtains 



dr 32 



i<J - 1) (7mm^CMB)""' X'^Q^^^iX, a) , 

(71) 



dt- 



e* - 1 



dt- 



tmin 



e* - 1 



duQ{X, t, u) 



duQiX,t,u), (72) 



where Q{X,t,u) is given by Eq. (j53p . and the boundary 
values are expressed by 



Un 



1 



mm "-mm 



1 



X 

^CMB7min , 



(73) 

(74) 
(75) 



and is given by Eq. ([54)1 . It is needless to men- 
tion that Qkn (^' ''') contains explicit 7min-dependences 
in contrast with Qkn {X, a) . 

In the Thomson limit x <^ 7inin, Eq. (j7ip is further 
simplified. Expanding Q{X, t, u) in powers of 1/it, the it- 
integration can be done analytically. Then, it is straight- 
forward to show that Eq. ([7T|) returns to Eq. ([TT|) as it 
should. 



B. Spectral intensity function for x < 7max region 

In the present subsection, we show the spectral inten- 
sity function with the Klein-Nishina formula for the case 
where the photon energy x is not far away from the upper 
boundary, namely x < 7max region under the following 
condition (denoted Case III hereafter): 

7min 'tiX < 7max ' ' ' CaSC III . (76) 

Under the condition of Case III, Eq. (|57)) becomes 
dl{uj) 



dr 



dsl{e ^xt) / d7Pe(7)-PKN(s,7) 



(77) 



Repeating the same procedure done in Sec. Ill, one 
finally obtains 

— -(Cr - 1) (7minC'CMBj X Qj^^(X,a), 



dr 32 



g^^(X,a)= / dt- 



(78) 



duQ{X,t,u), (79) 



where Q{X,t,u) is given by Eq. (j53p . and the boundary 
values are expressed by 



1 



4rniaxM: 



max tAmax 



^ max -J 

'Umax — ^ ^ J 

Tmax — ^'cMB7max , 



(80) 

(81) 
(82) 



and is given by Eq. (154^ . Again, it is needless to men- 
tion that (5|^[^(X, cr) contains explicit 7inax-dependences 
in contrast with QKNiX,a). 

Finally, for the extreme Klein-Nishina limit under the 
condition of 

X > 1, r 

max > 1 , (83) 
Eq. ([7^ is further simplified as follows: 

Q'^^i'^'^X, a) = 4C(2) [In X + C{X, a) + 1 - - Q] , 

(84) 
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where 
C{X,a) 



C7+1 



du 


1 


u 








1 - 


(/ r 




\ J- max / 



1 



In 1 - 



X 



n—1 ^ 

St 



1 + 2 In 2 
X 



(85) 



It is needless to mention that C{X, cr) returns to C{a) of 
Eq. deO]) for x < 7max (Case I). 



V. CONCLUDING REMARKS 

In NKl[l3, [IB], we have extended the formahsm ob- 
tained by Nozawa and Kohyama 16J to the case of high- 
energy electrons in the Thomson approximation. Assum- 
ing the power-law distribution for electrons, the analytic 
expression was derived for the spectral intensity function 
dI{uj)/dT. 

In the present paper, we have extend the formalism 
developed in NKI to higher photon energy (so called 
Klein- Nishina) region where the Thomson approximation 
breaks down. This work has been done on the basis of 
the recent work by Poutanen and Vurm[20j. where the 
general Klein-Nishina cross section has been calculated 
for the study of the inverse Compton scattering process. 

We first derived the analytic form for the redistribution 
function Pkn(s, 7). The obtained expression of Pkn(s, 7) 
is a natural extension of Pta(s,7) derived in the Thom- 
son approximation. Then, we have derived the expression 
for the spectral intensity function d/KN(w)/c?r, where 
the function is expressed by the double integral form. 
The form can be used in the entire Klein-Nishina en- 
ergy region (a; > 10^ eV). The function has an advan- 
tage that its computation can be done extremely fast 
compared with the full numerical calculation of multi- 
dimensional integrals which involve the numerical can- 
cellations among various terms. 

We have also derived approximate expressions for 
dlKNij^/dr, which agreed with Blunicnthal and 
Gould [7|. The expression for C{a) derived in the present 
paper has an advantage that it can be written in analytic 
forms for a = arbitrary fractional numbers. 

We have studied the knee structure of dlKN (i^) / dr . 
The physics origin on the knee structure is clear. The 
position of the knee is determined hy lu — to^c'*/A;bTcmb, 
which is the characteristic energy scale of the CMB pho- 
tons and electrons. The origin of the change of the slope 
in the spectrum is due to the change from the Thomson 
approximation to the extreme Klein-Nishina approxima- 
tion. It has been shown that the spectrum changes its 
slope at w « 10^^ eV from a = (tr — l)/2 to a = cr. The 
position of the knee is same as the well-known cosmic ray 



spectrum. However, the increment of the power index is 
Aa = 1.75 for a =2.5, which is quite large compared with 
Aa =0.4 of the cosmic ray spectrum. 

Finally, we have studied various cases of the boundary 
parameters (7min and 7max) of the electron distribution 
function, where the spectral intensity function depends 
explicitly on the 7niin and 7max parameters. For the ex- 
treme Klein-Nishina limit, we have shown an analytical 
expression for the spectral intensity function which takes 
into account for the finite size correction of the 7max pa- 
rameter. 
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Appendix A: Derivation of Q^^{X,a) 

We derive the expression of Q^^{X,(t) from Eq. ([5^ 
by taking the limit X ^ 1. According to Eq. ((55)) . one 
has It* ^ 1. The first term in Eq. ([53]) can be integrated 
analytically. One has 



du- 



1 



1 



(l + u)'^+2 a + 1 
The second term is 



(Al) 



du 



1 



u{l + u) 



cr+3 



du 



1 



(1 + uy+^ (1 + u)" 



-\niXt-J2-, 



(A2) 



i=l 



where n is a positive integer of the order of cr + 3. Equa- 
tion (|A2p can be shown as follows: First, define J„ by 



1 



du- , 

u(l + u)" 



(A3) 



(A4) 



into Eq. (jA3[) . and solving the recurrence relation for J„, 
one obtains 



Inserting the identity relation 
1 1 



m(1 + u)" u(1 + u)"-i (1 + m)" 



J n 

Ji 



n—1 ^ 

— ^ 7 



du— - 

u[l + u) 

In — KilnAXt. 



(A5) 



(A6) 
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One has 



du 



n—1 

ln4Xt-^-, (A7) 



Therefore, Eq. (jMl) is vahd. 

The third term can be neglected, and the fourth term 

is 



du 



1 



iXt. 



(A8) 



Inserting Eqs. (jATj) . (jA2|) and (jA7|) into Eq. ([52]), one 
obtains 

t 

Q'^^X, <j)^A dt- [hit + \nX + C{a)] , (A9) 

Jo e — i 



C{a) 



du 



{l + uY+^ (l + u)' 



9 ""^ 1 

^-l + 2hr2-y i 
tr + 1 i 

i—l 



(AlO) 



where we put u* — )■ 0. The ^-integration of Eq. (jA9P can 
be done analytically as follows: 



Qi=*(X,a)-4C(2)r(2) 



lnX + C(f7) 



C'(2) , r'(2) 



In deriving Eq. (|A11I) . we used 



dt , 





CG^)r(s): 



C(2) r(2). 

(All) 



(A12) 



dt- 







— -^[C(.)r(.)] 



= C(.s)r(.) 



Cis) n-s) 



(A13) 



Inserting C'(2)/C(2) = -6/^2 ^^^3 In fc/fc' = -Q, 
r'(2)/r(2) = I-C'e, and r(2) = 1 into Eq. (|ATT|1 . one 
finally obtains 



Q^'i^lX, a) = 4C(2) [InX + + 1 - - C,] 



(A14) 



where Ce = 0.5772 is the Euler's constant, and Ci 
0.5700. 



Appendix B: Analytic forms for C{a) 

We derive analytic forms of C{a) for specific cr- values. 
For example, the analytic form was derived for a — inte- 
ger cases in Eq. (pT|) . 



In order to derive the analytic expression for a = half- 
integer cases, let us introduce 



J, 



(i) - 



1 



du- 

u, u{l + u) 



n+i 



(Bl) 



where n > is an arbitrary integer, and ~ 1 /AXt <C 1. 
Rewriting Eq. (|B1I) with a new variable v — 1/(1 + 
one obtains the following recurrence relation: 



7(5) 7(5) 

'^TL •J-n-l 



2n- 1 

Equation (jB2p can be solved, and one has 



t(^) _ 7(3) 



n 



(I) 



In — = In l&Xt . 



(B2) 

(B3) 
(B4) 



Therefore, one obtains 



/ du r=lnl6Xt-2V . (B5) 

Ju, u(l + <+3 tt^*-! 



Thus, we have 



(l+1i)-+3 (l+u)"+5 

2 " 1 

l + 41n2-2V 

(T + 1 -^2^-1 



(B6) 



Finally, the analytic form for a = half-integers is obtained 
by inserting n=[a]-\-2> into Eq. (|B6I) as follows: 

2 1 
C(a)^-^-l + 41n2-2^^, (B7) 

i— 1 

where [cr] denotes the integer part of cr. 

Similarly, one can derive the analytic forms for a ~ q/p 
(irreducible fractional numbers). Let us introduce 



du 

IL, u{\ + u) 



(B8) 



where n is a natural number, and p and q are positive 
coprime integers. Rewriting Eq. (jBSP with a new variable 
V — 1/(1 -I- u)^/^, and solving the recurrence relation, one 
obtains 



1 



1-vP' 



(B9) 
(BIO) 



where v.^, ^ 1 ~ u^,/p. The elementary integral of 
Eq. (jBlOp can be done analytically. 
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We summarize the analytic expressions, for example, (vi) for cr = n + | 
in the fohowing cases: 
(i) for a = n 

2 '^1+^1 C(a) = — -4 l + 61n2 + - 

C7(a) = ^-l + 21n2-5:l, (Bll) H + f + l 2 



= 1 * H+3 



r-w ^1 -^-^^ln(^/2 + l)-8^ -1-, (B16) 

(n) tor 0- = n + g \/2 ~^ m — 6 

M+3 (vu)fora = n+| 

+ ^ + Sz^' ^^^^^ 

1=1 

2 TT 

/...x r , 1 C(a) = — o H-51n2 

(ill) for (7 = n + 4 ^ ^ + ^ + l 2 



t]+3 



(iv) for <T = n + 



4 

. , 4i-3 

(vmj lor a — n + 



C(cr) = — — ?5 l + 61n2-- ^ 2 1 , «, o 

M + i + 1 2 C(a) = ^-^^^-l + 61n2-- 

+^-V21n(V2+l)-8g^, (B14) - + ^2 ln(x/2 + 1) - 8 ^ , (B18) 

* i=l 

(v) for cr = n + ^ 

2 ['^1+3 where [a] denotes the integer part of a. Finally, their 

C(a) = — 1 + 4 In 2 — 2 > , (B15) numerical values are summarized in Table I. 

' [a] + i + l h^'-^ 
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TABLE I: Numerical values of C{a) ior a = n + a. 







n=0 


n=l 


n=2 


n=3 


n=4 


n=5 


n=6 


a 


=0 


0.8863 


-0.4470 


-1.030 


-1.397 


-1.664 


-1.873 


-2.046 


a= 


=1/8 


0.6159 


-0.5407 


-1.084 


-1.435 


-1.693 


-1.897 


-2.065 


a= 


=1/4 


0.3921 


-0.6267 


-1.136 


-1.471 


-1.720 


-1.919 


-2.085 


Q= 


=3/8 


0.2026 


-0.7061 


-1.184 


-1.506 


-1.748 


-1.942 


-2.104 


a= 


=1/2 


0.0393 


-0.7798 


-1.231 


-1.539 


-1.774 


-1.963 


-2.122 


a= 


=5/8 


-0.1038 


-0.8485 


-1.275 


-1.572 


-1.800 


-1.985 


-2.140 


a= 


=3/4 


-0.2306 


-0.9129 


-1.317 


-1.604 


-1.825 


-2.005 


-2.158 


a= 


=7/8 


-0.3443 


-0.9733 


-1.358 


-1.634 


-1.849 


-2.026 


-2.176 



